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Combining density functional theory (DFT) and embedded dynamical mean-field theory (DMFT) 
methods, we study the metal-insulator transition in R 2 R 2 O 7 (R= Y, Eu, Sm, Nd, Pr, and Bi) and 
the topological nature of the insulating compounds. Accurate free energies evaluated using the 
charge self-consistent DFT+DMFT method reveal that the metal-insulator transition occurs for 
an A-cation radius between that of Nd and Pr, in agreement with experiments. The all-in-all-out 
magnetic phase, which is stable in the Nd compound but not the Pr one, gives rise to a small Ir 4+ 
magnetic moment of ~ 0.5 hb and opens a sizable correlated gap. We demonstrate that within this 
state-of-the-art theoretical method, the insulating bulk pyrochlore iridates are topologically trivial. 


Interest in 5 d compounds has blossomed in recent 
years, as they provide a promising playground for in¬ 
teresting new physics^ - — arising from the interplay of 
the atomic spin-orbit coupling (SOC), which scales as 
Z 4 with the atomic number Z, and the electronic cor¬ 
relations, which are reduced due to the spatially more 
extended wavefunctions of the 5 d ions. Rare-earth py¬ 
rochlore iridates (R2R2O7) have drawn intensive at¬ 
tention in recent years^ because of their geometrically 
frustrated lattice, which favors the spin-liquid phased 
and the possibility to host various nontrivial topologi¬ 
cal phases in the bulker— thin films j 13 ^ 14 and domain 
wallsd Depending on the radius of the A-site cation, 
rare-earth pyrochlore iridates undergo a metal-insulator 
transition (MIT)jl& concomitant with a transition to an 
all-in-all-out (AIAO) magnetic stated - — Nevertheless, 
when the A-site cation is Pr or the more covalent Bi, the 
corresponding pyrocholore iridates are metallic down to 
very low temperature with no long-range magnetic order¬ 
ing fUfi 

Recently, inspired by pioneering work based on a 
generic tight-binding model^ there have been many stud¬ 
ies focusing on nontrivial topological phases in bulk py¬ 
rochlore iridates^ - — Based on calculations using the 
local density approximation (LDA) with Hubbard U 
(LDA+U) including SOC, Wan et al. predicted that 
Y2R2O7 can host nontrivial Weyl semimetal and axion 
insulator phases.— Using a simplified single-band Hub¬ 
bard model on the pyrochlore lattice, Go et al. con¬ 
ducted cluster-DMFT calculations^ and confirmed the 
LDA+U+SOC phase diagram, but with a magnetic con¬ 
figuration having the Ir magnetic moments rotated by 
90 ° from the AIAO ground-state configuration. 

Iridates are known to be beyond the range of ap¬ 
plicability of LDA+U calculations, since dynamical lo¬ 
cal correlations are crucial to explain their compli¬ 
cated electronic structure. Instead, the LDA+DMFT 
methocd correctly gives the boundary of the MIT in 
the Ruddlesden-Popper (RP) series, and the fine details 
of the ARPES measurements of S^IrO^d can only be 
understood by comparing with the LDA+DMFT spec¬ 
tral functions, which inlcude essentially exact correlation 
effects local to a given Ir atomd Moreover, the LDA+U 
method overestimates the stability of the insulating mag¬ 


netically ordered states, and the value of U, which is to 
some extent an adjustable parameter, often needs to be 
adjusted for each member of the series. Thus, a con¬ 
sistent treatment of the electronic correlations together 
with spin-orbit coupling at the LDA+DMFT level is nec¬ 
essary to shed light on the MIT in the pyrochlore iridates. 

In this work, we have carried out all-electron charge 
self-consistent LDA+DMFT— calculations. We have 
studied a series of pyrochlore iridates in both the param¬ 
agnetic (PM) and AIAO magnetic states, focusing on the 
metal-insulator transition with respect to the A-cation 
radius and the topological nature of the insulating states. 
We demonstrate that the MIT occurs in those compounds 
which can host an AIAO magnetic ground state. Be¬ 
cause of the large degree of geometric frustration in the 
pyrochlore lattice, a significant quasiparticle mass in the 
PM state is needed to destabilize the Fermi liquid forma¬ 
tion at the expense of magnetic long-range order. The 
tuning across the boundary between the Kondo-screened 
Fermi-liquid solution and the magnetic AIAO solution is 
achived with A-cation substitution. The most important 
consequence of the latter is the change of oxygen coordi¬ 
nates, which results in a slightly different Ir hybridization 
and effective Ir-^p bandwidth. Based on calculations of 
the wavefunction parities in a many-body context, we 
conclude that the insulating pyrochlore iridates are likely 
to remain topologically trivial. 

Our LDA+DMFT calculations were done using the 
projection/embedding implementation^ which avoids 
downfolding or model building, and adds dynamic local 
correlations to a set of localized quasi-atomic orbitals. 
The DFT part is implemented using the WIEN 2 k pack¬ 
aged All parameters of the calculation are identical to 
those used in our earlier work on RP iridates d We also 
performed LDA+U calculations using the VASP code 2 £ 
for comparison. Details are provided in the Supplemental 
Materials. 

For the series of R2R2O7 (R= Y, Eu, Sm, Nd, Pr, and 
Bi) compounds, we considered the PM and AIAO mag¬ 
netic solutions. We found that the magnetic AIAO solu¬ 
tion could be stabilized (that is, we could find a stable 
or metastable solution) only in the Y, Eu, Sm, and Nd 
compounds, and was always insulating (see spectral func¬ 
tions in the Supplement), while the PM solution could be 
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stabilized for all materials and was always metallic, con¬ 
sistent with experiments ] 16 i 2Q Naively, the MIT in the 
six R2IT2O7 compounds can be understood in terms of 
the bandwidth of the Ir -t2 g bands, as shown as inset of 
Fig .[lja). For instance, at the LDA level there is a critical 
value of about 2.6 eV for the Ir+23 bandwidth separat¬ 
ing the compounds with insulating and metallic ground 
states. A detailed analysis reveals that the bandwidth is 
closely correlated with the hybridization function, and its 
variation can be attributed to the change of the oxygen 
coordinates (see the Supplement). 

To determine the relative stability of the two phases, 
we computed the free energy differences between PM and 
AIAO states following Ref. [ 26 |, yielding the results dis¬ 
played in Fig. DJa). For comparison, the difference of 
LDA+U total energies is also shown, but reduced by 
a factor of 20 to fit in the plot. For each method, we 
used the same V as determined in our previous study 
of RP iridates^i The LDA+U energy difference is about 
400 meV on average, much too large to be compatible 
with the known ordering temperatures. This is hardly 
surprising since the local moments are zero in the PM 
state within LDA+U, which cannot describe fluctuating 
local moments. As a result, the energy of the PM state 
is severely overestimated. 

Experimentally, P^I^CU has a non-magnetic metallic 
ground stated correctly described by LDA+DMFT. In 
LDA+U by contrast, only BUL^Or is metallic, while the 
ground state of P^I^CU is an AIAO magnetic insulator. 
We also studied other possible magnetic phases within 
LDA+DMFT. We found that the three-in-one-out con¬ 
figuration can also be stabilized for some compounds, but 
is higher in energy than the AIAO state (e.g., by 23 meV 
for Y2R2O7). The two-in-two-out state, on the other 
hand, is not stable in this method. Therefore, the AIAO 
magnetic state is indeed the ground state of insulating 
pyrochlore iridates within LDA+DMFT. 

Our self-consistent LDA+DMFT calculations also give 
reasonable estimates of the free-energy differences be¬ 
tween PM and AIAO states as shown in Fig. H{a). The 
average free-energy difference is about 20 meV, which is 
about twice the typical magnetic ordering temperature 
(about 100 K) observed experimentally^ Moreover, the 
trend of the LDA+DMFT free-energy difference with 
respect to the A-cation radius follows loosely the ex¬ 
perimental magnetic ordering temperature^^ We note 
in passing that the difference in electronic entropy S 
between the two phases is very small. For example, 
in Nd2h*207 we obtained AS = 0.033 ks per Ir atom, 
which is similar to the experimentally determined AS = 
0.028 (see the Supplement). 

To shed more light on the origin of the MIT, Fig. DJb) 
shows the DMFT-predicted quasiparticle renormaliza¬ 
tion amplitude Z a = (1 — dlmE aa /dujn) -1 |o; n =(n where 
a labels a correlated orbital and £ indicates the self¬ 
energy on the Matsubara-frequency (uj n ) axis. The less 
active J e f? = 3/2 states are only weakly renormalized. 
Turning to the more active J e g = 1/2 states, these are 



FIG. 1: (Color online) Variation with respect to the A-cation 
radius of computed properties of R 2 R 2 O 7 (R= Y, Eu, Sm, Nd, 
Pr, and Bi). (a) Difference of free energies F between PM 
and AIAO magnetic states for R 2 R 2 O 7 . Solid circles (hollow 
squares) denote the values obtained using the LDA+DMFT 
(LDA+U) method, plotted as zero when the AIAO state can¬ 
not be stabilized. Inset displays the LDA band-width of the 
t 2 g bands; shaded region highlights the insulating compounds, 
(b) Quasiparticle renormalization amplitude Z a (see text for 
the definition) for the J e ff = 1/2 and J e ft = 3/2 states. Solid 
(dashed) lines denote the values in the PM (AIAO) states, 
with Z for the J e s = 1/2 orbital in the PM state of Sr 2 lr 04 
(red triangle) shown for comparison, (c) Ordered moment of 
Ir atoms. All LDA+DMFT calculations are done at 50K. 


renormalized by a factor of almost ten in the PM state 
for compounds having an insulating ground state. By 
contrast, the J e ff =1/2 quasiparticles in Bi2L*207 re¬ 
main quite light, consistent with recent ARPES measure¬ 
ment. 2 - Surprisingly, the J e ff— 1/2 orbital in Pr2lr207 is 
quite heavy (about five times the band mass), indicating 
that the correlations are considerably stronger than in 
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the marginal Mott insulator Sr 2 Ir 04 , whose is plotted 
as the red triangle in Fig. [TJb) . In this sense, the frus¬ 
tration in the pyrochlore lattice plays an essential role in 
shifting the MIT boundary towards stronger correlations 
compared to unfrustrated RP iridates. The frustration 
thus promotes metallicity by penalizing competing long- 
range magnetic order, an effect which is well captured by 
DMFT but less adequately by LDA+U, explaining why 
LDA+U predicts the Pr compound to be insulating. 

In short, the boundary of the MIT across the R 2 R 2 O 7 
compounds is determined by a competition between the 
formation of a quasiparticle band via the Kondo effect in 
the PM state, and the tendency toward long-range mag¬ 
netic order. Both can lead to an effective reduction of the 
free energy, and are favored by enhanced quantum fluctu¬ 
ations. In general, the occurrence of a MIT depends cru¬ 
cially on many factors, including dimensionality and frus¬ 
tration. For example, in Sr 2 Ir 04 the two-dimensionality 
reduces the Neel ordering temperature, but the short- 
range order above the Neel temperature still preserves 
the charge gap in the excitation spectrum. Frustration 
also reduces the tendency to long-range order, but it pro¬ 
motes metallicity and allows very narrow quasiparticles 
to be observed in the PM state before the long-range or¬ 
der opens the charge gap. This is why Pr 2 Ir 2 C >7 remains 
in the paramagnetic bad-metal phase even though the 
correlation strength is much larger than in Sr 2 Ir 04 . 

For correlated Fermi liquids, the quantum fluctuations 
give rise to fluctuating moments in the PM state, which 
statically order when the magnetic energy gain is suffi¬ 
cient to overcome the gain available from band forma¬ 
tion. According to Fig. DJc), the magnitude of the or¬ 
dered magnetic moment for Y 2 Ir 2 07 is about 0.57 /ib- 
This is a bit larger than the experimental upper bound 
of 0.5 fiB obtained by neutron scattering^ which can be 
attributed to the fact that the spatial fluctuations, which 
are not fully accounted for in LDA+DMFT, are likely 
to further reduce the ordered magnetic moment. In the 
AIAO insulating state, there are no states at the Fermi 
energy, but the slope of the self-energy at zero frequency 
still gives some measure of the importance of quantum 
fluctuations. As shown by the dashed lines in Fig. [TJb), 
the quantum renormalization effects are greatly reduced 
in the AIAO state as compared to the PM state. 

We observed that, like the effective J e ff = 1/2 states 
in the RP iridates the J e ff = 1/2 wavefunction in the 
pyrochlore iridates also strongly deviates from the rota- 
tionally invariant SU( 2 ) point. This is mostly due to 
the trigonal crystal field induced by compression of the 
Ir06 octahedra. The trigonal crystal-field splitting in 
our calculations is derived from a very large energy win¬ 
dow, and is on average about 0.39 eV (see the Supple¬ 
ment), in good agreement with resonant X-ray scattering 
measurements^ Correspondingly, the resulting orbital- 
to-spin moment ratio is about 1.3 (see the Supplement), 
showing a significant deviation from the SU{ 2 ) limit of 
two *21 Nevertheless, the effective J e ff = l/2 and J e ff = 3/2 
bands are still well separated in energy, which facili- 


TABLE I: Parity analysis of eight eigenstates around E f 
in the PM and AIAO magnetic states, corresponding to 
Figs. [2{a-b) respectively. ‘MULT’ denotes the multiplicity 
in accounting for the eight time-reversal invariant momenta. 
The parities are given in order of increasing energy eigenval¬ 
ues, with vertical bars indicating Ef . 
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tates our analysis of the topological character. Note that 
the values of local U and J parameters on Ir atoms in 
our LDA+DMFT calculations are the same for the py¬ 
rochlore and RP iridates *21 confirming the more universal 
nature of local Coulomb repulsion across similar materi¬ 
als when the screening by hybridization effects at high 
energy are allowed in the calculation. 

We now turn to the topological properties of the in¬ 
sulating pyrochlores. Fig. [2fa) shows a fictitious non¬ 
magnetic band structure of Nd 2 Ir 2 07 , together with 
the irreducible representations for the J e ff = 1/2 states 
around Ep at the eight time-reversal-invariant momenta 
(TRIM). This band structure is obtained by neglecting 
the dynamic part of the LDA+DMFT self-energy, but us¬ 
ing the self-consistent LDA+DMFT charge density, with 
time-reversal (TR) symmetry imposed. We note that 
these bands are quite similar to the LDA bands. 

Regarding the irreducible representations, the space 
group of the pyrochlore iridates is nonsymmorphic. Thus, 
depending on the choice of inversion center, the four en¬ 
ergetically equivalent L points are separated into one L' 
and three L points, with opposite parities for the L' and 
L points as shown in Fig. [2fa). Moreover, at each of 
the three X points all states are four-fold degenerate by 
symmetry, with each degenerate group comprising two 
even-parity and two odd-party states. Finally, all eight 
J e ff = 1/2 bands at the T point have even parity. The 
situation is summarized in the PM column of Table HI It 
follows that if a global band gap could be opened so as 
to extend the shaded region of Fig. [2^a) throughout the 
Brillouin zone, corresponding to a half-filled J e ff = 1/2 
manifold, seven out of 16 of the occupied Kramers pairs 
at the eight TRIM would be odd-parity ones. Since this 
number is odd, it would generate a strong-topological- 
insulator phased It has been suggested this might be 
achieved for Pr 2 Ir 2 07 , for example, by applying a strain 
along the (111) direction . 31 If the TR symmetry would be 
weakly broken in this phase, the system would become an 
axion rather than a strong topological insulator, but we 
shall use the term “topological insulator” to cover both 
cases. 

When the DMFT self-energy is taken into account, the 
TR symmetry is spontaneously broken and a global gap 
opens up. However, the topological character may be dif¬ 
ferent from that described above, because of the drastic 
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Self energy (%) 

FIG. 2: (Color online) (a) The LDA-like band structure 
of Nd 2 lr 2 07 with irreducible representations indicated at 
the time-reversal invariant momenta; numbers in parenthe¬ 
sis denote degeneracies, (b) Color map: spectral function of 
Nd 2 lr 207 in the AIAO magnetic state. Dotted lines: band 
structure of the effective topological Hamiltonian (see text for 
details). The parities of eight J e ff= 1/2-derived bands are in¬ 
dicated by “+” or ” at the L point, (c) Evolution of the 
eigenvalues of the effective topological Hamiltonian at L with 
mixing fraction of the DMFT self-energy. Inset zooms in on 
the critical region indicated by the red square. 


mixing of bands that results. Following the theory of in¬ 
teracting topological phases ; 32 ! 33 the topological indices 
can be obtained by inspecting an effective single-particle 
Hamiltonian defined as H e ^ = Hq + E(u; = 0), where Hq is 
the TR-invariant Bloch Hamiltonian and E(u; = 0) is the 
DMFT self-energy at zero frequency (i.e., at the Fermi 


energy), which carries all the TR-breaking terms. The 
resulting band structure is shown in Fig. [2jb), together 
with the full LDA+DMFT spectral function. Since TR 
symmetry is broken, the total parity of all occupied eigen¬ 
states at the TRIM should now be counted^ the system 
will be a trivial or axion insulator if the total number is 
4 n or 4n + 2 respectively, with n an integer. 

The parities of the J e s = 1 / 2 -derived eigenstates at the 
eight TRIM are also shown for the AIAO state in Table HJ 
At the r point all states remain even-parity even with the 
addition of the self-energy. At the X points, each four¬ 
fold degenerate state splits into two doublets, and we find 
(see the Supplement) that including Y(cj = 0 ) does not 
induce any band inversion or otherwise change the order¬ 
ing of parities, so that the number of occupied odd-parity 
states remains equal to two at each X point. The most 
significant change occurs at the L (Z/) points, where the 
TR-symmetry breaking splits each doubly-degenerate L 4 
(L 4 ) state into two singlets. Some of these cross the gap 
as the self-energy is turned on, as shown in Fig.[2jc), with 
the result that there is an exchange of parity between oc¬ 
cupied and unoccupied states as shown in Table HJ As a 
result, there are a total of twelve odd-parity occupied 
states at the TRIM. Similar behavior is also observed in 
the other insulating R 2 IT 2 O 7 compounds studied. Thus, 
after properly including the DMFT self-energy, we find 
that the insulating R 2 I 12 O 7 compounds are topologically 
trivial insulators. This is in agreement with previous 
LDA+U result^ in the large -U limit. 

Nevertheless, in Refs. 0 and [Hit was argued that non¬ 
trivial axion insulator and Weyl semimetal phases occur 
with smaller values of on-site U. This is not supported 
by our DMFT calculations. For instance, for Y 2 R 2 O 7 , 
which is the most insulating of the six compounds, the 
PM solution is more stable than the AIAO magnetic state 
when the on-site U value is reduced to 4.0 eV (see the 
Supplement). Moreover, as shown in Fig. [2fc), we ob¬ 
serve that the parity exchange between occupied and un¬ 
occupied states at the L and L' points happens before a 
global band gap is opened. Therefore, we think it is very 
unlikely that a topological insulating phase can be found 
in the bulk pyrochlore iridates. 

In summary, our calculations provide a clear picture of 
the origin of the MIT in pyrochlore iridates and its vari¬ 
ation with respect to the A-cation radius. Moreover, our 
parity analysis of the many-body effective J e ff = 1/2 wave 
functions reveals that when insulating magnetic phases 
appear in these pyrochlore iridates, they are very likely 
to remain topologically trivial. 
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Note: As we were finalizing this manuscript a preprint 
appeared^ which presents results for Y 2 lr 20 y that are 
similar to our own in some respects. However, we do 


not agree that the magnetic moments on the Ir atoms in 
Y 2 H 2 O 7 can be larger than 1.0 /jlb (see the Supplement 
for a detailed discussion). 
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